Spin Hall effect in two-dimensional p-type semiconductors in a magnetic field 
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We calculate the spin Hall conductivity driven by Rashba spin-orbit interaction in p-type two- 
dimensional semiconductors in the presence of a perpendicular magnetic field. For a highly confined 
quantum well, the system is described by a fc-cubic Rashba term for two-dimensional heavy holes. 
The eigenstates of the system can be described by Landau spinor states. First we consider the 
conventional spin Hall conductivity. The contribution of the interband transitions to the Kubo- 
Greenwood formula gives the density dependent intrinsic spin Hall conductivity, which approaches 
its universal value a% y — 9e/Sn for weak spin-orbit coupling and low Fermi energies, in agreement 
with previous work. However two intraband contribution terms cancel this effect leading to zero 
conventional spin Hall conductivity. Adding the torque dipole contribution to the definition of spin 
current, we also study the effective spin conductivity. This is shown to be proportional to the total 
magnetization plus surface terms which exactly cancel it for small spin-orbit coupling. If in low 
magnetic field the intraband transitions evolve to vertex corrections, the fact that both effective and 
conventional spin Hall conductivities vanish is unexpected. This suggests that the zero magnetic 
field limit of the model is singular. 

PACS numbers: 72.25.-b,72.25.Dc 
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I. INTRODUCTION 

The goal of spintronics is to inject and manipulate the 
spin of charge carriers in semiconductors. Spin Hall effect 
(SHE) has been proposed as an efficient source of spin 
current in semiconductors. Two different mechanisms of 
SHE are proposed^ In the extrinsic SHE the skew scat- 
tering of electrons (driven by an electric field) from im- 
purities with spin-orbit (SO) dependent scattering poten- 
tial leads to opposite spin accumulation in the transverse 
edges of the sampled On the other hand, the intrinsic 
SHE arises from the SO term in the single-particle Hamil- 
tonian of paramagnetic materials and is present even in 
the absence of impurities. Intrinsic SHE is predicted for 
two major systems: Murakami et al£ have predicted a 
dissipationless spin current in p-type semiconductors de- 
scribed by the Luttinger model 4 of spin-3/2 holes. Sinova 
et alA have calculated spin Hall conductivity (SHC) for 
a two-dimensional electron gas with k-linear Rashba 6 SO 
coupling. 

Two seminal experiments have reported the observa- 
tion of SHE. Kato et alJ have used Kerr rotation to 
observe spin accumulation in n-type GaAs. Wunderlich 
et al& have observed the effect in a clean 2D hole system 
using circularly polarized light. It is generally believed 9 - 
that in the former experiment the measured SHE is ex- 
trinsic while in the latter experiment (due to very low 
impurities concentration) it is an intrinsic effect piLlLiSi 

When impurities are present their contribu- 
tion to intrinsic SHE is crucial. Extensive 
theoretical, 13 i 14 i 15 i 16 i 17 i 18 i 19 and numerical studies, 20 ! 21 ! 22 
show that for the model of 2D electrons with fc-lincar 
Rashba SO interaction, the SHE does not survive the 
thermodynamic limit if disorder induced vertex correc- 
tions are included. This is because in this system the 



spin current is the time derivative of the magnetization 
and vanishes when the magnetization saturatesi 1 ^ 6 " On 
the other hand for the Luttinger model of a hole gas in 
which the SO coupling depends nonlinearly in momen- 
tum, the vertex correction is calculated to be zero within 
the ladder approximation and for short-range scattering 
potential of the impurities ^LiLiSi Adding Rashba SO 
interaction in highly confined quantum wells, makes the 
system be described by fc-cubic Rashba SO termsiS 3 . The 
spin Hall conductivity of this system in the clean limit 
approaches the universal value of a = 9e/8n for weak 
SO coupling^ 5 . The vertex corrections of this model are 
calculated jiSiii* 3 ^ and confirmed numerically ^StSLSL to 
be zero. 

Rashba has used an external magnetic field on a 2D 
electron system with fc-linear Rashba SO term to show 
that the zero spin current is an intrinsic property of 
the system and persists in the presence of impurities. 28 
This system is described by two series of Landau levels 
(spinors) separated by the SO dependent gap. The stan- 
dard Kubo- Greenwood formula has been used to calcu- 
late the SHC. It turns out that the interband transitions 
contribute to the universal intrinsic SHE while including 
the intraband transitions cancels the effect. In the zero 
magnetic field limit the latter contributions evolve into 



vertex corrections 
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The relation between the conventional definition of 
spin current (as the expectation value of the product of 
the spin operator and velocity) and spin accumulation is 
not obvious. In systems with SO interaction the spin cur- 
rent is not conserved. The spin density continuity equa- 
tion contains the torque contribution which adds a torque 
dipole term to the conventional spin current term^ 4 *^ 7 . 
The resulting effective spin current has been introduced 
as the proper definition of the spin current. The cor- 



2 



responding effective spin Hall conductivity is shown to 
have opposite sign with respect to the conventional SHC 
for fc-linear and fc-cubic Rashba models^ 

Parallel to recent work, 25 ' 2 £ we consider a 2D hole gas 
model in a perpendicular magnetic field in the absence 
of Zeeman coupling and find its exact eigenstates and 
eigenvalues. The Kubo-Greenwood formula is then used 
to calculate both the conventional and effective SHC of 
the model for different values of SO coupling and Fermi 
energies, including both interband and intraband contri- 
butions. We show that unexpectedly both effective and 
conventional SHC go to zero for small value of SO cou- 
pling and then take a finite value as the SO coupling 
increases. Although calculations of vertex corrections 
for different types of scatterers have actually appeared 
recently; 3 ^^ more detailed comparison and analysis will 
be desirable especially in the presence of magnetic field. 

II. THE MODEL 

In 2D hole systems, the Rashba spin-orbit splitting 
(due to the structure-inversion asymmetry) of heavy hole 
states is of third order in in-plane wave vector fc as com- 
pared to Rashba spin splitting of electron states which is 
linear in For a very confined quantum well and not 
too high densities, the Hamiltonian of heavy holes can be 
written as 

Ho = + + 7o^ 3 (fc^+ + k\a-\ (1) 
2m 2m 

where k± — k x ± \k yi er^ = (a x ± ia y )/2. The Pauli 
matrices a represent spin-3/2 heavy holes with effective 
mass m, and 70 is the Rashba SO coupling which without 
the lose of generality we assume to be positive^ Out of 
two energy branches 

e(fc) = h 2 k 2 /2m ± ~f Q h 3 k 3 , (2) 

the lower energy branch is bounded from below only for 

70 < l/(3mkh). (3) 

For low Fermi energies (or low hole densities no), this 
condition is satisfied and the SHC is calculated in the dc 
limit neglecting the impurity vertex corrections^ When 
the impurity scattering energy H/t is comparable with 
the SO splitting energy 7 (2mey) 3 / 2 , the SHC goes to 
zero. In the opposite limit of a clean system H/t <C 
7o(2me/) 3 / 2 <C £/, the dc SHC is found to be dependent 
on 7on D : 

Here k^(k^ ) is the Fermi momentum for the upper 
(lower) energy branch Recent theoreticalji2*ii and 
numerical r2i2SiSi work shows that the vertex corrections 



to the SHC of the above model are zero for short-range 
disorder interaction. This is in contrast to the fc-linear 
Rashba model in which the vertex corrections lead to 
complete cancellation of the spin current. 

We consider the above homogeneous model in a ho- 
mogeneous perpendicular magnetic field B = —Bz, in 
the absence of Zeeman coupling^ The kinetic momen- 
tum P = p + -A satisfies the commutation relation 
[P y ,P x ] = imhw Cl where u> c = eB/mc. We use the sym- 
metric gauge A = — ir x B. The Hamiltonian is 

H = ^ + ^+7o(^a++P 3 a-) 
2m 2m 

= hu c ^(a f a+ i) +7(a 3 cr+ +a^a-)^j (5) 

in which the bosonic operators a = P-/\/2mhuj c and 
= P + 1 \J2mhio c describe Landau levels. The dimen- 
sionless parameter 7 is defined by 7 = 'yo\/8m 3 fiiL> c . 

Similar to the fc-linear Rashba model in a perpendic- 
ular magnetic field ^^^^ the exact eigenvalues and 
eigenstates of the above model can be found. The SO 
term couples every state (</> n ,0) to the (0, ^,1+3) leading 
to two new A = ± spinor states: 

,/.A = ( cos 0^_a \ 
Vn { sin^ n+1 J ' 

E*/hu c = n + 6/ cos26*, (6) 

in which S = -3/2, tan20+ = ^/n{n 2 - 1) 7 /S, 6~ = 
®n ~ tt/2 where n > 2. The states <fioi,<fiil and 4>2l are 
not coupled to any other state, i.e for n < 2 we have 

e = k = 0. 

We restrict the upper limit of SO coupling to 7 < 
l/(3y / ri/) where ej — huj c rif is the Fermi energy. This 
condition is equivalent to equation |J3J. As long as this 
condition is satisfied, two A = + states do not cross each 
other. Energy levels are shown in Fig.^ we can see that 
A = — states never cross each other (for 7 > 0). 

III. THE EFFECTIVE AND CONVENTIONAL 
SPIN CURRENT 

Using the commutation relation of P x and P y the ve- 
locity can be written as 

v x = [x,H]/ih= [P y ,H]/(imfkj c ) 

v y = [y,H]/ih = -[P x ,H]/(imhu c ). (7) 

The same relation also holds for the linear SO sys- 
tem in magnetic field, 28 but the spin-dependent part 
of the velocity is different in these two models. While 
in the hole system velocity shows a d-wave symmetry, 
v x = P x /m + 3j(P x P v + P y P x )a x + 3 7 (P 2 - P%)&» , for 
the electron system we get v x = P x /m + ja y in which 
the 7-dependent part is a constant. The symmetry of 
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torque dipole density. The dc SHC is 

= l3 % 2 Y {(Xn\P x \X'n')(X'n'\P y a z \Xn) 

nn , XX 

(Xn\P x \X'n')(X'n'\[P y a z + a z P y , V}\Xn) ^ 

2{E*-E») 1 J 

in which the first term o^?* (respectively the second term 
<J xy ) 0n the right hand side is the contribution of the 
effective spin current (respectively torque dipole density) 
to the conventional SHC. 



FIG. 1: The energy levels E% of eq.@ as a function of 7 for 
3 < n < 20. Note that each energy state E£ is plotted as long 
as E% > E^_ 1 . The inset shows all states in that window, 
including A = — states that disperse upwards and so do not 
cross each other. 



the velocity is directly reflected in the vertex corrections. 
For the 2D hole system in the absence of magnetic field 
this leads to zero vertex corrections after integration over 
azimuthal angle <j) — ^W- ipy/Px)^^ 

The conventional definition of spin current for spin-3/2 
holes is 



-(v x a z +a z v x ) = 3P x a z /2m. 



(8) 



Note that in contrast to the £;-linear Rashba model, here 
the spin current operator can not be written as a com- 
mutator like J* ~ [H, a x ]m& This is the profound reason 
behind the vanishing of spin current in fc-linear Rashba 
model for any scalar random potential. 

Relations Q and (JSJ) show that 
(Xn\v y \X'n / )(X'n'\J z \Xn) is antisymmetric and so 
we can use the following definition of conventional 
SHC^ 
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(Xn\v y \X'n')(X'n'\J z \Xn) 



(E*-E*y-h 2 w 2 



(9) 

where the sum is restricted by > £/ and E*> < £/. 
The Landau degeneracy factor mu) c / {2irh) is also in- 
cluded. 

Using the relation J z = | — [Py<r z ,Ho] the conven- 
tional spin current can be written as 



V 



2mfku. 



HPyV Z 



-a z P y ,H] 



[P y a z + a z P yi V]) (10) 



where V is the SO term in (JSJ. The first term in this 
formula is a full time derivative and represents the effec- 
tive spin current of the model (jSJ) ■ By simple algebra the 
second term can be written as P y & z + cr z P y which is the 



IV. RESULTS AND DISCUSSION 

Relations Q and JSJ indicate that the transition from 
the state n is allowed only to the state n' — n± 1. In the 
energy spectrum ©, the A = + state has lower energy 
than A = — levels (because we have defined 8 < and 
7 > 0). For a given Fermi energy the first unoccupied 
A = + state n + is larger than the first empty A = — state 
n~ i.e always n + > rT +3. Therefore according to Fig. [3 
there are two contributions to the spin Hall conductivity: 
the interband transitions occur from any state if~ with 
71 > nT to the state V'nii with nil < n + , and the 
intraband transition from ib + > to ?b + , . and from ib~_ 

^n+— 1 r n 

to ^~__ r 

To calculate the effective SHC in 1)110 we apply the 
same method used by Rashba^ namely we write 



a z xy s (0) = 3i- 



7Tmco c h 2 

( V (-l,n|P x |A'n')(AV|P,a 2 |-l,n) 



n+-l 



n—n 



+(-l,n + \P x \ + l,n+ - l)(+l,n+ - l\P y cr z \ - l,n+) 
+(+l,n+|P x | + l,n+ - l)(+l,n+ - l\P y <r z \ + 
-<-l,n+ - 1|P X | - l,n + ){-l,n+\P y o z \ - l,n+) 
-<-l,n+ - 1|P B | + l,n+)(l,n+|P a a z | -l,n+ - 1)). 

(12) 

In deriving this equation we have used the fact that the 
transition ip~ — > Vv+i * s cance l e d by Vv+i ~~ * V'n ■ The 
equivalence of 1)1211 and 1)11) 1 can be seen from Fig. [21 
The sum of the last four 'surface' terms in 1)12) 1 renders 
a constant —3/2 for the fc-cubic Rashba model using JBJ. 
On the other hand taking into account the x—y symmetry 
of the model, the first term in 1)12)1 can be written in 

the symmetric form J2n=n-(~^' n \[Pxi P y ]c z \ — L™)/2- 
Finally we get the following relation for the effective 
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FIG. 2: (Color online) Left panel shows possible interband 
(dotted line) and intraband (dash-dotted line) transitions of 
the model. The solid (dashed) line shows filled (empty) states. 
For convention the A = — 1 states are plotted lower. The 
right panel shows the equivalence of equations 1121 and Hill . 
Two subtracted x transitions and two added • transitions 
are 'surface' terms in 1121 . 



SHC: 



3e 



^ S (°) = ^( E (-l,»k«|-l,n>/2 + 3/2) 



34r(3/2- V cos26>„/2). 



(13) 



The second term is the total magnetization of the system 
in a magnetic field. 

For a given physical system the SO parameter 70 
may depend on the intrinsic and applied electric fields^ 
The Fermi energy depends linearly on the hole density 
ri£). Moreover "i^/nj = (2m) 3 ' 2 7oy / £/ 1S independent of 
the applied magnetic field. In the accumulation layer 
of GaAs — Alo.sGao.sAs one obtains jy/nj « 0.3 for 

rip = 10 10 cto~ 2 while in a 200A-wide quantum well with 

,-2 
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n d = 10 cm we get "fy/nj ~ 10" 

We study the SHC as function of this parameter. 
When plotted against this parameter all graphs of SHC 
for different Fermi energies (hole densities) will be scaled 
to one graph. To begin with, the value of 7 is lim- 
ited by 1 < 7^^ 2 < rif which corresponds to h/r c < 

7o(2me/) 3 / 2 < £/, i.e the same conditions as in the free 
hole system in which the relaxation time r is replaced 
by the cyclotron time r c = l/u> c . For this value of 7 the 
mixing angle (9 A approaches its maximum value — 7r/4 and 
the conductivity (fill) can be approximated as: 



-9^(l/7Vn+-l/ 7 Vn- + l). (14) 



In Fig. 01 we have plotted the effective SHC (jl 1|> as a 
function of jy/nj. For larger value of jy/nj the effective 



SHC <r%.'y is negative and finite. The negative sign of the 



effective SHC for larger jy/nj is the result of the original 
choice for positive conventional SHC. 29 

For lower values of "fy/nj both the effective and con- 
ventional SHC approach zero. In Fig. 0] we have plotted 
the effective SHC for lower values of Jy/nj . Going to 
the zero limit the typical Shubnikov-de Haas oscillations 
appear in the conductivity. Notice that for smaller n/ 
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FIG. 3: (Color online) The effective spin Hall conductivity 
as a function of 7y / nJ = 70 y/ej for different values of Fermi 
energy. All graphs scale to one if plotted against yy/nj. For 
weak SO coupling or low Fermi energies the conductivity ap- 
proaches zero but it shows Shubnikov-de Haas oscillations. 
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FIG. 4: (Color online) The effective spin Hall conductivity 
for very low value of 7y^nJ ~ lOy/^J- The higher the nf val- 
ues the smaller is the amplitude of oscillations, as one would 
expect. The inset is exactly the same graph when plotted 
against -yuf 3 ^ 2 which shows that the highest jump in SHC is 
^ a xy = 9e/(87r) and it happens at the same value of 772/ 3//2 
for any n/. 



the amplitude of the oscillations is larger, as one would 
expect for the quantum limit. 

It is important to note that the same behavior can be 
observed for effective SHC in the linear Rashba model in 
a magnetic field. Actually the formula (|11|) is valid for 
/c-linear Rashba model (up to the factor of 9 which comes 
from spin-3/2 holes and fc 3 momentum dependence). In 
Fig. [S] we compare the effective SHC of these two mod- 
els. There are two differences between the two systems. 
Firstly in the fc-linear Rashba model with Vso = c«(fc x c) 
the effective SHC remains zero for larger value of SO 
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FIG. 5: (Color online) The distinct intra and interband tran- 
sition contributions to the effective SHC for very low value of 
7y / nJ in the fc-cubic model (left panel). The right panel is 
the same for the fc-linear model in which two contributions 
are zero independently. 
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FIG. 6: (Color online) The interband contribution to the 
conventional SHC as a function of 7-ynJ = 70-y/eJ for the 
fc-cubic model. For weak SO coupling or low Fermi energies 
the conductivity approaches its universal value a% y = 9e/87r. 
When the intraband contribution is included, the total con- 
ventional SHC is zero for small 7y / fTJ and goes to negative 
values as "iJnf increases. 



coupling a while in the fc-cubic model it goes to finite 
negative values (see Fig. [2J ■ However the deviation from 
zero in the latter case is because we approach the limit of 
unbounded negative energies while for the linear model 
the typical value of SO coupling is very far from the up- 
per limit a > y/nj of the unphysical region. Second and 
more importantly, in the linear model both the interband 
and intraband contributions to the effective SHC are zero 
independently while for the fc-cubic model these two con- 
tributions are finite but opposite in sign and cancel each 
other out. 

For comparison we also calculate the conventional SHC 
@ . A straightforward replacement of the wave functions 
and energies in @ renders 



(|15|l can be approximated by 



'(0) 




(16) 



which using hwcU^ 1 = f?k\j2vn renders exactly Q). 
Here, as 7 decreases the conventional SHC approaches 
its universal value u% y — 9e/8ir^ Now we include also 
the intrabandr transitions, consisting of the two terms 
near the Fermi surface. Their contribution is negative 
and cancels the conventional SHC for the lower values of 
J-JnJ, as it is shown in the inset of Fig. H3 



n, XX' 



(n — 1) COS 



2 aX' 



J n+1 ' 



_,2 aX 



^-(n + 2)sin 2 ^; i 



.2 a x 



1 + 5/ cos 29% +1 - 5/ cos 20* 
(n - 2) cos 2 0*'_ x cos 2 0* - (n + 1) sin 2 0*'_ x sin 2 0* 



1 + 5/ cos 20* - 5/ cos 20*'_ 1 



(15) 



in which the first (respectively second) term describes 
the transition n — > n + 1 (respectively n — > n — 1). We 
first include only the interband contribution to the con- 
ventional SHE to compare it to the system of free holes 
in the absence of magnetic field. The result is plotted in 
Fig. and it is actually the same as the SHC for a free 
2D electron gas^ For the range of "f^/rTJ chosen here, 



V. CONCLUSION 

We have studied the 2D hole gas system with fc-cubic 
Rashba SO interaction in a perpendicular magnetic field 
in the absence of Zeeman splitting. The exact eigen- 
states and eigenvalues are derived and used in a Kubo 
formula to calculate the effective and conventional SHC 
of the model. For the conventional SHC the interband 
transitions reproduce the intrinsic spin Hall conductivity 
which approaches the universal value of a^ y — 9e /8n for 
weak SO coupling. 25 However when the negative intra- 
band contributions are included the conventional SHC 
becomes negative and goes to zero for weak SO coupling. 

We also consider the effective spin current in which the 
torque dipole is included. We show that the correspond- 
ing effective SHC shows Shubnikov-de Haas oscillations 
around zero for small SO coupling, if both interband and 
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intraband contributions are included. Higher n/ results 
in smaller oscillation amplitude. For stronger SO cou- 
pling the effective SHC becomes negative. 

Comparing the fc-cubic model with the fc-linear model, 
the main difference is that in the former the intraband 
and interband contributions to the effective SHC are fi- 
nite but opposite in sign and cancel each other. In the 
latter each contribution is zero independently. 

The fact that for low values of SO coupling both the 
effective and conventional SHC go to zero is unexpected. 
For the fc-linear model it is well-known that the intra- 
band transitions evolve into the vertex correction which 
is proved to be zero. For the fc-cubic model, although 
it is well knowniSiii*iS that the vertex corrections are 
zero in the case of short-range s-wave scatterers, the re- 



lation between the intraband transitions and the vertex 
corrections is not straightforward. For the effective SHC 
ca. (|13|l suggests that the zero magnetic field limit is sin- 
gular in this model. However more detailed calculations 
of the vertex and higher order corrections may help to 
ascertain this behavior. 

The above results would be modified for an inhomo- 
geneous system and for the ac SHC a* y (uj), because the 
intraband contributions are less stable and are reduced 
by frequency and scattering^ One expects that the in- 
traband contributions will diverge at u) = oj c after which 
they will vanish for uj c <C u). 
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